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The full three-dimensional dispersion of the 7r-bands, Fermi velocities and effective masses are 
measured with angle-resolved photoemission spectroscopy and compared to first-principles calcula- 
tions. The band structure by density-functional theory underestimates the slope of the bands and the 
trigonal warping effect. Including electron-electron correlation on the level of the GW approxima- 
tion, however, yields remarkable improvement in the vicinity of the Fermi level. This demonstrates 
the breakdown of the independent electron picture in semi-metallic graphite and points towards a 
pronounced role of electron correlation for the interpretation of transport experiments and double- 
resonant Raman scattering for a wide range of carbon based materials. 



Recently graphene has been investigated as a proto- 
type system to address basic questions of quantum me- 
chanics (relativistic Dirac fermions) as well as for 
high speed semi-metal field effect transistors in emerg- 
ing nanoelectronic devices [4l. Many of these results 
are based on its peculiar electronic properties, i.e. an 
isotropic and linear dispersion close to the Fermi level 
{Ep). In low dimensional and strongly anisotropic sys- 
tems correlation effects play a crucial role in understand- 
ing and describing the electronic band structure. Kinks 
in the quasiparticle (QP) dispersions and lifetimes were 
observed and interpreted as band renormalization due 
to electron-phonon [5 and electron-plasmon |^ interac- 
tions and band structure effects [7] . Its electronic proper- 
ties are also very sensitive to stacking and the number of 
layers [5]. In bilayer graphene a gap that could be tuned 
by doping was observed [9J. For few layered graphene, 
the parent compound, graphite, is the key to understand- 
ing these new phenomena. Interlayer coupling in an AB 
stacking sequence leads to the formation of electron and 
hole pockets responsible for the semi-metallic character 
in graphite. The linear dispersion is broken and only if 
we have an AA stacking the linear dispersion remains. 
Nevertheless, at the H point of graphite [2] the 7r-band 
dispersion is close to linear and has been interpreted as 
Dirac- Fermion like. Much less is known about the quan- 
titative description of electron-electron correlations in 
these graphitic systems. Angle-resolved photoemission 
(ARPES) studies indicated that local density approxi- 
mation (LDA) gives a dispersion that is too flat and a 
scaling has to be applied in order to fit the experimen- 
tal dispersion of few-layer graphene and graphite. For 



the scaling, values of -10% ^ E] and -20% HH [T3] 
have been reported in the literature. Correlation effects 
can cause the recently discovered quantitative shortcom- 
ings [14] of the double-resonant model [15] for Raman 
scattering in graphene and graphite. Furthermore, it is 
important to know the exact dispersion, because it 
is responsible for the conductivity perpendicular to the 
graphene layers. 

In this letter we report on a combined ARPES and 
theoretical ab-initio QP study of the three dimensional 
TT band structure and the Fermi surface in graphite sin- 
gle crystals. ARPES is best for studying correlations 
since it probes the complex self-energy function which 
contains the electronic interactions. We elucidate the 
full electronic QP dispersion perpendicular to the layers 
and show directly the importance of the effects of elec- 
tronic correlation and the influence on transport and res- 
onance Raman in these layered materials. Experiments 
were done at BESSY II using the UE112-PGM2 beamline 
and a Scienta SES 100 analyzer yielding a total energy 
resolution of 15 meV and a momentum resolution better 
than 0.01 A"^. Natural graphite single crystals with AB 
stacking were cleaved in-situ to give mirror-like surfaces 
and were measured within 12 h after cleavage in a vacuum 
better than 10^^*^ mbar. The samples were mounted on a 
three axis manipulator that was cooled down by liquid He 
to 25 K. The three-step model for ARPES is employed 
for data analysis [16|. After a dipole transition from the 
valence band to an unoccupied intermediate state, the 
electron travels to the sample surface where it is scat- 
tered into a free electron state outside the sample. When 
the electron leaves the sample, its wavevector component 
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FIG. 1: (color online) (a) Raw EDC data for the cut through H. Cuts at special ky (values in A~^) are indicated by red lines. 
Photoemission intensity of cuts through the 3D BZ of graphite at (b) H, (c) kz = 0.22A~^ and (d) K point. Along with the 
photoemission intensity, we show results of LDA (□, red line) and GW (A, blue line) calculations. The inset in (b) shows a 
magnification of the TB-GW calculation for the small hole pocket and the gap A. The photon energies used were (b) 100 eV, 
(c) 25 eV and (d) 83 eV. The intensity in (c) is multiplied by a factor of 2 in order to visualize the low intensity lower tt band. 



in the direction perpendicular to the sample surface is not 
conserved which makes 3D band mapping a challenging 
task. We assume that the intermediate state dispersion 
is a parabola which is shifted downwards in energy with 
respect to the vacuum level by the inner potential. We 
determine the inner potential by finding the maximum 
splitting between the two valence tt bands, which occurs 
at kz = 0, i.e. K. This procedure yields an inner poten- 
tial of Vb = 16.4 ± 0.1 eV and allows one to accurately 
determine the full dispersion by changing the photon 
energy |16j . Deviations from a parabolic dispersion of the 
intermediate state would change the kz assignment as de- 
scribed in [17l UHl [19] . However, as shown in detail below 
and in agreement with previous experiments [21 [12] for 
the electronic states of graphite close to the Fermi level 
this error is small. 

The calculations of the electronic dispersion are per- 
formed on two levels. First, we calculate the Kohn-Sham 
band-structure within the LDA to density-functional the- 
ory (DFT) [50]. Wave-functions are expanded in plane 
waves with an energy cutoff at 25 Ha. Core electrons 
are accounted for by Trouiller-Martins pseudopoten- 
tials. In the second step, we use the GW-approximation 
[2T1 [25] to calculate the self-energy corrections to 
the LDA dispersion [53]. For the calculation of the di- 
electric function e(w, g) we use a 15x15x5 Monkhorst- 
Pack k sampling of the first BZ, and conduction band 
states with energies up to 100 eV above the valence band 
(80 bands). The calculations are converged to 10 meV. 
Since the GW calculations are computationally expen- 
sive a tight-binding fit to GW (TB-GW) was performed 
in order to rapidly calculate the QP bands. 

We now turn to the comparision of the ARPES data 



to the calculated band structures. In Fig. [T] we show 
raw EDC data and cuts of the tt bands that intersect 
the corners of the BZ at k^ = 0.46 A~^ (very close to H 
point), at kz = 0.22 A^^ and at fc^ = (K point). The 
cuts are done along the ky axis [see coordinate system 
in Fig. [2j^a)]. It is clear that for the cut through the H 
point, the bands are almost linear and only one tt valence 
band can be seen, whereas at the K point, the degener- 
acy is lifted and the bands are parabolic. The energy 
separation between the two valence tt bands is getting 
smaller when moving away from K in kz direction. The 
LDA and GW calculation data along the measured cut 
are also shown in Fig. [l] LDA underestimates the slope 
of the bands while the self-energy corrections from the 
GW calculation yield a better agreement. In the inset of 
Fig. [T]jb) we show the gap A at H obtained by GW. It 
is clear that the dispersion is not linear close to H and 
the concept of Dirac Fermions breaks down. This is im- 
portant for transport properties, however for estimation 
of the band maxima a linear extrapolation above Ep is a 
good approximation due to the extremely high curvature 
at the maximum of the hole like band at H. 

We can obtain the Fermi surface of the hole like 
band with a band maximum above Ep by the following 
method: we fit a linear dispersion to the EDC maxima 
in a region below Ep and then extrapolate it to find the 
crossing of the extrapolated band with the Fermi level. A 
dispersion very close to linear is predicted by calculation 
[inset of Fig. [l|^b)] as the functional form of the band at 
H. From the crossing of the linear dispersion with Ep, 
we can then obtain an estimate for the cross section of 
the hole pocket at H. Such a procedure is needed in or- 
der to determine the maximum of the sparsely populated 
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FIG. 2: (color online) (a) Graphite BZ with the high symme- 
try points, reciprocal lattice vector c and the coordinate sys- 
tem. The red (green) curves around K{H) are the envelopes 
of the measured electron(hole) pockets, (b) depicts the inter- 
polation of the 3D Fermi surface to experimental data, (c) 
Obtaining band maxima: crosses denote symmetrized EDC 
data and the green line an extrapolation. The hole pocket 
cross section with Ef is denoted by green circles. The arrow 
denotes the band maximum at H. (d) photoemission inten- 
sity along KH direction for photon energies from 22 eV to 
105 eV and GW calculations (blue line). The five colored 
points along KH denote for which we evaluate the EDC 
maxima, (e) Raw EDCs for the denoted by colored points, 
(f) The experimental EDC maxima and the calculated LDA 
and GW dispersion. 

band and has been used in the literature 2 . For points 
with kz < 0.5c, we choose a quadratic fit as the func- 
tional form around the band maximum (c is the length 
of the reciprocal lattice vector in direction) . Assuming 
electron-hole symmetry for the sparsely occupied band, 
we can estimate the cross sections with Ep. From this 
data and some higher energy points the Fermi surface 
could be interpolated with the standard tight-binding 
Hamiltonian pS]. 

The pockets we obtained are shown in Fig. |2|a) along 
with a sketch of the BZ and the coordinate system and 
reciprocal lattice vector c. A detailed magnification of 
the tiny electron and hole pockets is shown in Fig. [2|b) 
with the aspect ratio enlarged by a factor 30 for eas- 
ier viewing. The surfaces around K and H points are 
the electron and hole pockets, respectively. The vol- 
umes inside the Fermi surface based on experiment yield 
carrier densities of Ue = 2.3 x 10^^ cm~'^ electrons and 
Uh = 1.8 X 10^^ cm~'^ holes. In Fig.[2]jc) we illustrate the 
procedure to obtain the tt band maxima at H which is lo- 
cated above Ep. In Fig. [2Fd) we show the photoemission 
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FIG. 3: (color online) Equi-energy contours of the photoe- 
mission intensity (logarithmic scale) around the KH axis for 
fcz = 0.22 corresponding to Fig. 1(c). Experimental 

EDC maxima are denoted by crosses. Lines depict TB-GW 
(blue, A) and LDA (red, □) calculations, respectively. Trig- 
onal warping and the strong asymmetry in the intensity can 
clearly be seen. 



intensity along the KH direction for photon energies be- 
tween 22 eV and 105 eV. In Fig. ^e) the raw EDC data 
for selected cuts along KH are depicted and in Fig. 
the EDC band maxima (folded back to the first BZ) are 
compared to theory. The experimental results and theory 
have the same trend, which is that the lower tt valence 
band comes close to the upper one when moving from K 
to H point. The experimental splitting between the tt 
valence bands is generally better described by the GW 
approximation than by LDA. 

In Fig. [3] the trigonal warping effect is investigated. In 
order to accurately compare to calculations we evaluate 
the maxima of the EDCs for equi-energy contours in the 
kj:,ky plane. We have chosen fc^ = 0.22 A~^, which lies 
approximately halfway between K and H. The experi- 
mental and calculated GW and LDA contours for three 
energies 0.1 eV - 0.3 eV are shown in (a)-(c). LDA not 
only underestimates the slope of the bands but also a 
simple scaling does not work because a different scaling 
would be needed for F and M directions. Although the 
GW approximation does not describe the experimental 
data in all directions, it is always closer to the experiment 
than the LDA. Interestingly, it is in perfect agreement 
with the experiment at H point while the agreement is 
worse when going away from H. This discrepancy might 
be due to the GW approximation and resolved by inclu- 
sion of higher orders in the calculation of the self-energy 
correction in a future work. 

A strong asymmetry around the KH axis in energy 
(trigonal warping effect) and photoemission intensity is 
observed. The trigonal warping effect of the ARPES data 
is underestimated by the LDA and again better described 
in the GW approximation. Compared to graphene we 
find that trigonal warping is higher by a factor of ~ 1.5 
when considering the distances in the two high-symmetry 
directions for an equi-energy contour of 0.1 eV in Fig. [3] 
The strong asymmetry in the photoemission intensity 
seen in this figure is attributed to the dipole matrix ele- 
ment for the transition from a valence band state to an 
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Method 


VF (HA) 


VF (HL) 


energy 

separation (eV) 


A (eV) 


ARPES 


1.06 ±0.1 


1.05 ± 0.1 


0.71 ± 0.015 




GW 


0.99 


1.02 


0.71 


0.037 


LDA 


0.77 


0.87 


0.62 


0.022 



study (mg is the free electron mass) [23. For the scat- 
tering time r we use a value of t ~ 200 fs from pump- 
probe experiments [28 . We get a DC conductivity a — 



^r/ml + nhC^T/ml of a = 3.9 x 10^ 



The 



TABLE I: Fermi velocities at H (in 10® m/s), valence band 
energy separation at K and A obtained by experiment and 
from calculations. 



unoccupied state 25 eV above Ep 

We now turn to a quantitative assessment of electron 
correlation effects and their implications for Raman and 
transport studies of graphite and related materials. The 
experimental and calculated values for the energy sepa- 
ration of valence bands at K, the gap A at H and the 
Fermi velocity vp are listed in Table |l] We obtain the 
Fermi velocities of the linear band in HA and HL di- 
rection in a region up to 1.5 eV below Ep. Since the 
trigonal warping effect has a minimum in this kz plane, 
the differences in these two directions are small. The 
LDA values for vp are far off and the differences of 
in HA and HL directions (trigonal warping effect) are 
strongly overestimated. On the other hand the GW is in 
excellent agreement with the experiment. 

The dispersion of D and G' double resonant Raman 
(DRR) bands [TS] with laser energy is proportional to 
vphlvp, where vph is the slope of the transversal opti- 
cal phonon branch going through K (H). The electronic 
dispersion in kz broadens the DRR dispersion and with 
the present data DRR spectra for 3D graphite can be 
calculated taking into account the full electronic disper- 
sion. For a given phonon dispersion we thus expect the 
present data to reduce the slope of the calculated DRR 
bands since the vp from our experiment are blue-shifted 
with respect to the LDA calculation. 

In addition, the gap A shows a large difference in val- 
ues with a GW value that is five times the TB parameter 
that was originally fit to a magnetoreflectance experi- 
ment |25j . Such a deviation can be attributed to stack- 
ing faults such as AA stacking, which has A = 0. For 
a sample consisting of AB stacking with many stacking 
faults the effectively measured A would be lower explain- 
ing the experimental value of A = 0.008 eV in pyrolytic 
graphite [25]. A finite value of A causes a breakdown of 
the linear bands at the H point. Their shape becomes 
parabolic with a very large curvature and thus a small 
absolute value of the effective mass. 

From the 3D Fermi surface the transport proper- 
ties of graphite can be calculated from a Drude model 
where the carrier densities and the effective masses can 
be directly taken from the ARPES experiment. For 
the kz averaged electron and hole masses, we obtain 
TO* = O.IOtoq and m^^ = — 0.04too, respectively that is 
in good agreement with a previous magnetoreflectance 



values agree nicely with transport measurements that 
yield 2.5 x 10^ n'^cm.-^ ^ and 2.5 x 10^ n-^cm-^ [25] 
when considering that the literature values differ by an 
order of magnitude. 

In conclusion we have performed ARPES of graphite 
and compared the measured QP dispersion to ab-initio 
calculations. We have found that the band dispersions 
are better described by the GW approximation, however 
around the K point that agreement becomes poorer and 
many-body approximation schemes beyond GW may im- 
prove this result. Moreover, for this system the concept 
of Dirac Fermions breaks down. This highlights the im- 
portance of electron correlations which renormalize the 
electronic dispersion resulting in a blue shift of the gap at 
H and an increased vp. A direct comparison of a scaled 
LDA QP dispersion to transport and Raman experiments 
is not possible due to an anisotropy of the scaling needed 
for the LDA around the KH axis. TB-GW allowed us 
rapid calculation of the QP bands. 
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